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ABSTRACT 


In  order  to  determine  the  Fourier  transform  of  a  qua si- per iodic  time 
series  (linear  problem),  or  the  power  spectrum  of  a  stationary  random 
time  series  (quadratic  problem),  it  is  desirable  that  data  be  recorded 
without  interruption  over  a  long  time  interval.  In  practice,  this  may 
not  be  possible.  The  effect  of  regular  interruption  such  as  the  day / 
night  cycle  is  well  known.  We  here  investigate  the  effect  of  irregular 
interruption  of  data  collection  (the  ^breaking**  of  the  window  function) 
with  the  simplifying  assumption  that  there  is  a  uniform  probability  ^ 
that  each  interval  of  length  ,  of  the  total  interval  of  length  T  «  Nt, 

yields  no  data.y-^ 

For  the  linear. case  we  find  that  the  noise-to-slgnal  ratio  will 

—1  —2 

have  a  (one-sigma)  value  less  than  e  if  N  exceeds  p  (l-p)e 
For  the  quadratic  case,  the  same  requirement  is  met  by  the  less  restric¬ 
tive  requirement  that  N  exceed  p  *(l-p)e~*. 


It  appears  that,  if  four  observatories  spaced  around  the  earth  were 
to  operate  for  25  days,  each  for  six  hours  a  day  (N  "  100),  and  if  the 
probability  of  cloud  cover  at  any  site  on  any  day  is  20%  (p  ■  0.8), 
the  r.m.s.  noise-to-signal  ratio  is  0.25%  for  frequencies  displaced 
from  a  sharp  strong  signal  by  15  mHz.  The  noise-to-signal  ratio  drops 
off  rapidly  if  the  frequency  offset  exceeds  15  mHz. 


EXAMINATION  OF  TIME  SERIES  THROUGH  RANDOMLY  BROKEN  WINDOWS 

I.  INTRODUCTION 

In  many  astrophyslcal  problems  one  Is  concerned  with  the  study  of 
time  series.  It  often  happens  that  the  property  of  particular  interest 
is  the  spectrum  of  the  time  series.  In  principle,  one  may  determine  a 
time  series  to  a  prescribed  accuracy  by  making  measurements,  without 
interruption,  over  a  sufficiently  long  time  interval.  In  practice,  the 
length  of  time  over  which  the  variables  may  be  measured  will  be  limited. 
Moreover,  measurements  may  necessarily  be  interrupted  (or  otherwise  impaired) 
for  one  reason  or  another.  The  relationship  of  the  spectrum  determined 
by  limited,  interrupted  measurements  to  the  intrinsic  spectrum  has  been 
the  subject  of  many  investigations,  as  recently  reviewed  by  Deeming  (1975). 

If  the  original  time  series  is  denoted  by  x(t),  one  may  regard  the 
measurements  y(t)  as  being  determined  by 

y(t)  »  f (t)  x(t) ,  (1.1) 


where  f(t)  is  the  "window  function."  We  regard  x,  y  and  f  as  being  simple 
scalar  functions  but  the  procedure  may  be  generalized  to  replace  x,  y  by 
vectors  and  f  by  a  tensor. 

We  use  the  Fourier  transform  notation 


x(t) 

*(«) 


(1.2) 

(1.3) 


where  the  limits  of  integration  are  to  be  taken  to  be  - 00  to  +  «  if  other 
limits  are  not  explicitly  specified. 


If  we  are  interested  in  determining  x(co),  the  Fourier  transform  of 


the  time  series  x(t),  then  we  may  use  the  relation 


y(o>) 


j. 


do)'  f  (CD1)  x(0)  -  0)') 


(1.4) 


to  relate  the  Fourier  transform  of  the  measured  time  series  y(t)  to  that 
of  the  original  time  series  x(t). 

We  are  interested  in  the  possibility  that  f(t)  may  be  regarded  as  a 
random  variable,  expressible  as 


f(t)  *  F(t;a^,ci2 . ajj)  » 


(1.5) 


where  a^»a2*  *  * ,aN  are  independent  random  variables  with  specified  expec¬ 
tation  distributions.  By  the  central  limit  theorem  (Papoulis,  1965),  we 
expect  that  the  random  variable  f  (or  its  Fourier  transform  f)  will  have 
a  distribution  approximately  Caussian  in  form  if  N  is  not  a  small  number, 
so  that  an  adequate  representation  of  f  would  be  given  by  its  mean  <  f  > 
and  its  standard  deviation  0(f). 

If,  on  the  other  hand,  x(t)  is  a  random  time  series,  we  will  be 
concerned  with  the  autocorrelation  function  Rx(t),  defined  by 

Rx(t)  -  <x(t')  x(t’+t)>,  (1.6) 


and  the  power  spectrum  of  the  time  series,  defined  as  the  Fourier  transform 

of  R  (t): 
x 


R  (t) 
x 


(1.7) 


Sx(o» 


(1.8) 
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On  noting  that 


<x(2)x(oj')>  ■  Sx(<o)  6(<u  +  o»')  (1.9) 

and  evaluating  < y (<u)  y  (oi')>  ,  we  may  verify  that 

V“>  -  Ji*'  "<»’)  S,  <«-«•>  (1.10) 

where 

w(o»)  -  f  (cj)  f  (-0J).  (1.11) 

Clearly  the  function  w(u)  represents  the  capability  of  the  measurement 
process,  described  by  the  "window  function"  f(t),  to  determine  the  power 
spectrum  S^(oi) .  The  function  w(cu)  may  be  expressed  in  terms  of  the  inde¬ 
pendent  random  variables 

w(<u)  ■  W((OJ  oi^Q^, .  •  • ,  o^)  •  (1.12) 

Once  again,  unless  N  is  a  small  number,  we  expect  that  the  distribution 
of  w  will  be  approximately  Gaussian  so  that  it  may  be  characterized  by 
its  mean  value  <w>  and  standard  deviation  o(w). 

This  article  was  prompted  by  a  problem  related  to  the  determination 
of  normal  modes  of  oscillation  of  the  sun,  as  determined  by  measurement  of 
the  photospheric  velocity  field.  Measurements  have  been  presented  by 
Deubner  (1975)  and  by  Rhodes  et  al.  (1977),  and  their  theoretical  inter¬ 
pretation  discussed  by  Ulrich  and  Rhodes  (1977)  and  by  Ulrich  et  al.  (1978). 
For  optimum  determination  of  the  power  spectrum  of  the  velocity  field 
(expressed  as  a  function  of  wave  number),  it  is  clearly  desirable  to  make 
observations  without  interpretation  over  as  long  an  interval  as  possible. 
Away  from  polar  regions,  observations  from  a  single  station  are  interrupted 
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by  Che  day-night  cycle  which  leads  to  unacceptable  aliasing  of  the  data. 
Observations  made  from  a  spacecraft  in  polar  orbit  would  obviously  yield 
un-aliased  data  of  higher  quality  and  higher  frequency  resolution.  Obser¬ 
vations  made  from  the  south  pole  during  austral  midsummer  can  lead  to 
several  days  of  uninterrupted  observation  and  to  still  longer  intervals 
with  occasional,  irregular  interruption.  It  is  also  possible  to  select 
three  or  four  stations  around  the  earth  which,  in  the  absence  of  any  cloud 
cover,  could  give  continual  coverage  of  the  sun  for  many  weeks.  However, 
one  must  anticipate  that  some  of  the  data  would  be  lost  by  cloud  cover. 

It  is  clearly  desirable  that  one  should  be  able  to  make  some  esti¬ 
mate  of  the  accuracy  with  which  oscillation  modes  may  be  determined  when 
it  appears  possible  to  observe  the  sun  over  a  long  Interval  of  time  losing 
some  blocks  of  time  because  of  cloud  cover.  The  purpose  of  this  article 
is  to  develop  a  model  which  enables  us  to  address  problems  of  this  type. 
After  presenting  a  few  general  formulas,  we  shall  simplify  the  problem 
considerably  by  supposing  that  observations  are  made  over  a  large  number 
N  of  equal  time  Intervals,  each  of  length  T  ,  so  that  the  total  time 
interval  T  is  given  by 


T  -  NT  .  (1.13) 

With  certain  additional  simplifying  assumptions,  we  shall  consider  the 
statistical  properties  of  the  functions  f(co)  and  w(o>)  which  are  repre¬ 
sentative  of  "randomly  broken"  window  functions. 
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II.  MATHEMATICAL  MODEL 

In  the  case  that  the  window  function  f(t)  is  expressible  in  the  form 
(1.5),  in  terms  of  a  number  of  random  variables,  we  wish  to  study  the  distri¬ 
bution  of  the  functions  f(oj),  w(a>),  entering  equations  (1.4)  and  (1.10). 

We  suppose  that  the  distribution  of  the  variables  to  is  given  by 

the  probability  function  P(a^ . a^)  such  that  P(a^,  ...»  ct ^)dot^  ...,  da^ 

is  the  probability  of  finding  a^  in  the  range  to  a^  +  da^,  etc.  Then  the 
expectation  value  of  the  quantity  f(<u)  is  given  by 

<  f  (cu)>  -  JdNa  PN(a)  F (cojo^ . aN),  (2.1) 

where  dNa  denotes  da.  ...,  da„,  and  P_.(a)  denotes  P(a-.,  ...,  a„) .  If  we 

INN  x  N 

use  the  following  notation  for  the  variance  of  a  complex  variable  of  a 
complex  variable  z, 

a2(z)  -  02(zr)  +  a2(z±),  (2.2) 

where  z^  and  z^  are  the  real  and  imaginary  parts  of  z,  then  noting  that 
f  (-  cd)  Is  the  complex  conjugate  of  f(co),  we  see  that 


a  (f (co))  ■  <f(oj)  f  (-o>)>  -<f(<o)><f(-co)>  .  (2.3) 

The  first  term  on  the  right-hand  side  may  be  evaluated  from 

<w(oj)>  =<f(co)  f  (-oj)>  -  ^jNot  PN(a)F(a);a1,...,aN)  F  (-«;  c^,...,^)  (2.4) 

We  see  that  equation  (2.4)  also  gives  the  expectation  value  of  the 
"window  spectrum"  w(u)  which  appears  in  equation  (1.10)  and  is  appropriate 
for  the  discussion  of  stationary  random  time  series.  The  variance  of  this 
function  is  given  by 


a2 (w)  ■  <w(a>)  w  (- ai)>  -  <w(co)>  <w(-a>)> 


(2.5) 
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where 


Hence 


<\an>‘  P2  +  <P-P2>5»„.  «-13> 

where  <5^  is  the  Kronecker  function.  In  evaluating  the  quantity  given 

by  (2.6),  we  need  to  evaluate  the  expectation  value  of  a  a  a  a  .  By  con- 

m  n  p  q  J 

sidering  the  various  possibilities  (m,n,p,q  all  different:  two  of  them 
the  same,  etc.)  we  find  that 


<a  a  a  a  > 

mnp  q  ^ 


4  3  4 

p  +  (p  -  p  )  (6 


mn 


6  + 

mp 


mq 


+  6  + 
np 


nq 


6  ) 
pq 


where  6 

mnp 


+  (p2  -  3p3  +  2p4) (6 _ 

npq 

+  (p  -  4p^  +  6p3  -  3p4) 
1  if  *  •  n  ■  p  otherwise  0, 


+  <S  +  5  +6  ) 

mpq  mnq  mnp 


mnpq 


(2.14) 


and  is  defined  similarly. 
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III.  EVALUATION  OF  MODEL 


For  simple  (non-random)  time  series,  equation  (1.4)  gives  the  relation¬ 
ship  between  the  Fourier  transforms  of  the  original  and  measured  time  series. 
In  this  context,  the  properties  of  the  random  window  function  f(t)  may  be 
characterized  by  ^f(co)^  and  cr2(f). 

On  substituting  the  form  (2.8)  into  (1.5),  we  find  that 


F(o;<x^,  . . . ,  oc^) 


TO)  Lj. 

n  ■  1 


a 

n 


(3.1) 


On  using  (2.1)  and  (2.12),  we  obtain 

1  1  |ia)(t  +t  ) 

<f(o>)>  *  T  p  sine (y  0)  T)  e  1 

where  sine  0*0^  sin  0 . 

On  using  (2.4),  we  find  that 

<w(co)>-  (i)2  T2  |p2sinc2  (io>T)  +  N_1p(l-  p)sinc2( ycox) 

Hence,  using  (2.3),  we  obtain 


(3.2) 


(3.3) 


a2(f)  -  2  T2  N"1  p(l  -  p)  sine2  avr)  .  (3.4) 

For  evaluating  the  effects  of  "breaking"  of  the  window  function,  it 
is  convenient  to  normalize  the  standard  deviation  with  respect  to  the  maxi¬ 
mum  value  of  f((o),  which  is  the  value  at  w  =  0.  Accordingly,  we  introduce 
the  definition 


I1(0)) 


0  (f  (0))) 

<mr>  • 


For  the  case  under  consideration,  this  has  the  form 

I1(co)  -  N-1/2p'1/2  (l-p)1/2  |sinc|o)T|  . 


(3.5) 


(3.6) 
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For  discussion  of  the  properties  of  randomly  broken  windows  in  the 
study  of  stationary  random  time  series,  it  is  necessary  to  evaluate  the 
mean  value  and  standard  deviation  of  w(co) .  The  former  is  given  by  equa¬ 
tion  (3.3).  The  first  term  inside  the  brackets  has  the  same  form  as  arises 
in  the  non-random  case  (p  =  1).  The  second  term  represents  a  change  in  the 
mean  spectrum,  so  it  is  convenient  to  introduce  the  symbol  for  the 

ratio  of  the  additional  term  to  the  maximum  value  of  the  principal  term: 

A2  =  N-1  p-1  (1-p)  sine2  (icuT)  .  (3.7) 

On  writing  equation  (2.6)  in  the  simpler  form 

<w(co)w (—  oo) >  =  <F(co;a)F(co;ci)F(-  to;  a)F(-  co;  a)>  (3.8) 


and  using  equations  (1.3)  and  (2.8),  we  see  that 

<w(co)w(-co)>  =  / — “r,  <a  a  a  a  >(« 

li~..}4_W.  m  n  p  q  \ 


(zttco/  mnpq 

/  -icot 
x  ^e  1 


icot  iwt 
le  m  -  e 


/  icot  icot 

m-11  /  n  n-11 

\  “e  J 


-icot 


-  e 


■4 


-icot 


P -Ll\e  *  -  e 


-icot 


-4 


(3.9) 


On  using  equation  (2.14),  we  see  that  this  may  be  expressed  in  the  form 

<w(co)w(-co)> - — ,  !p\  +  (p3-p4)(E  +E*+  4E  ) 

(2ttco)^  11 


+  2(p2  -  3p3  +  2p4)(E4  +  E*) 


(3.10) 


+  (p  -  4p2  +  6p3  -  3p4)  E5j 


where 
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<w(a>)  w  <-<o)>  -  jp  sine  [j  wT| 


8 incut 
s  incur 


+  2(p3-p4)n“2  sinc2(-|  ut|  ainc2|-^  urj 
+  ^ (p3-  P4)  N-1  ainc2(-|ut)alnc2(-|  ur| 


+  4(p2-  3p3  +  2p4)  N-2  ainc2|io)Tj  ainc2|^toT^ 
+  |p  -  4p2  +  6p3  -  3p4)  N*3  8inc4|j  on^  |  . 


On  using  equations  (2.3)  and  the  definition 


ZjiO))  - 


a  (w(u)) 

( W  ' 


we  find  that 


2(w)  ■  |2  N  ^p  ^(1-p)  sine2  sine2  ^urj 
+  2N  2p  ^(1-p)  sinc2^  utj  sinc2^ urj 


sinaiT 

sinaix 


N  2p  2(3  -  lOp  +  7p2)  sinc2|y  ut|  sinc2^y  urj 

| 1/2 

+  N  3p  3(1  -  4p  +  6p2  -  3p3)  sinc4(-r- cot)(  . 
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IV.  DISCUSSION 


We  see  from  the  preceding  section  that  the  effect  of  a  random  "break¬ 
ing"  of  the  window  function  is  to  produce  an  aliasing  of  any  signal.  This 
effect  is  described  by  the  function  Z^  or  by  and  Z^  as  given  by  equa¬ 
tions  (3.6),  (3.7)  and  (3.15). 

For  the  "linear"  problem,  as  described  by  equation  (1.4)  the  mean 
Fourier  transform  of  the  window  function,  as  given  by  equation  (3.2)» 
has  the  same  form  as  it  does  in  the  non-random  case,  although  it  is 
reduced  by  a  factor  p.  The  standard  deviation  is  characterized  by  Z^, 
defined  by  equation  (3.5)  and  given  by  equation  (3.6). 

It  is  convenient  to  introduce  the  notation 


--1 


-1 


2T  ,  a)T  =  2t  . 


(4.1) 


-1 


On  noting  that  sine  0  <  1  for  9^1,  and  that  |sinc  9|  <  0  for  0^1, 


we  see  that  Z^  <  S^  where 


S1(ai)  -  N~1/2p"1/2U-P)1/2,  0)^cot 
S1(o>)  -  N"1/2p_1/2(l-p)1/2(<ox/to)  ,  0)  ^o>t. 


(4.2) 


Hence  the  aliasing  is  most  severe  within  the  range  of  a  few  times  oi^  of 
the  strongest  signal. 

We  may  infer  from  the  above  restriction  the  minimum  number  of 
intervals  necessary  to  ensure  that  Z^  is  below  an  assigned  level  e  for  a 
given  value  of  p. 

We  see  from  (4.2)  that  we  require  N  >  Nowhere 

Nx  -  (l-p)p"1e"2.  (4.3) 
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If,  for  instance,  p  ■  0.8  and  we  require  that  E^ <  0.05,  N  must  be  at 
least  100. 

For  the  "quadratic"  problem  in  which  we  are  determining  the  spectrum 
of  a  stationary  random  time  series,  the  aliasing  is  described  by  the 
functions  A2(co)  and  E2(<o)  given  by  equations  (3.7)  and  (3.15). 

We  see  from  equation  (3.3)  that  the  second  term  in  brackets  is  smaller 
than  the  envelope  of  the  first  term,  and  so  may  be  neglected,  for  co  <  co^  , 
where 

a>A  -  N1/2(l-p)1/V1/\  (4.4) 


For  co  >  ci>A ,  the  second  term  produces  a  "tail"  to  the  principal  contribu¬ 
tion  to  <w(o))>  .  By  an  argument  similar  to  that  leading  to  equation  (A. 2), 
we  find  that  &2  <  D2>  where 

D2(co)  -  N"1p‘1(l-p)  ,  co<  coT,  ) 


(4.5) 


D2(w) 


n‘1p_1(1-p)  (a>T/0))2 


co  >  01 


Now  consider  the  four  terms  in  the  second  bracket  in  equation  (3.15). 

It  is  clear  that  the  third  term  may  be  ignored  by  comparison  with  the  first 

since  it  has  a  similar  dependence  on  co  but  includes  an  extra  factor  N 

The  second  term  may  also  be  ignored  by  comparison  with  the  first;  the 

extra  factor  N  ^  (sin  coT)/(sin  cor)  has  a  maximum  value  of  unity,  and  an 

-1/2 

RMS  value  of  order  of  N 

In  comparing  the  fourth  term  by  comparison  to  the  first,  we  see  that 
the  ratio  is  given  by 


R 


sinc2|yCoT| 

sinc2(-jC0tj 


(4.6) 
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where 


Q(P) 


1  -  3  p  +  3p' 

,  2 

2p 


(4.7) 


It  is  easily  verified  that  Q(p)  <0.5  in  the  range  0.5  <  p  <  1.  Hence 
R  has  a  maximum  value  of  order  0.5  at  values  of  co  for  which  cot  ■  2nir.  We 
find  that,  when  R  is  averaged  over  frequency,  it  varies  with  N  as  N  *. 

Hence  we  may,  to  sufficient  approximation  ignore  the  fourth  term  and  so 
replace  (3.15)  by 


£  2(co)  * 
We  find  that  T 2  < 
S2(co) 

S2(oj) 

s2(w) 


2l/2N-l/2p-l/2(1  _  p)1/2|8tnc|Ia)Tj||slnc|i.,JJX)|. 

S2 ,  where 


.  i1  Wlftd  -Dl/2. 

-  21/2N~1/2p_1/2(l  -  p)1/2 

-21/2N-3/V1/2(l-p)1/2 


w  <  a)T  » 


coT  <  CO  <  co^  , 


(O  >  <0^  . 


(4.8) 


(4.9) 


We  see  that,  for  the  same  values  of  N  and  p,  the  maximum  value  of  is 
1/2 

2  times  larger  than  the  maximum  value  of  E^ .  Hence  the  minimum  number 
of  Intervals  necessary  to  ensure  that  is  below  an  assigned  level, 
for  various  values  of  p,  is  twice  the  corresponding  value  of  N^,  given  by 
equation  (4.3). 

However,  the  quadratic  case  is  more  complicated  than  the  linear  case 
in  that  Z2(co)  is  more  complicated  than  E^(co),  and  Aj(co)  is  nonzero  (whereas 
A^(co)  is  zero  and  has  been  neglected). 


On  noting  that  the  dominant  term  of  (3.3)  (that  which  survives  in  the 

2 

nonrandom  case  that  p  *  1)  varies  as  (u^/oj)  for  o>  >  a>T>  we  see  that  I2(o>) 
is  less  than  the  tail  of  the  dominant  term,  and  so  may  be  neglected,  for 

oj  <  ojj.  ,  wht:. 


(0 


Z 


2-l/2N1/2p1/2 


(1-p) 


(A. 10) 


On  the  other  hand,  we  find  that  S2(o>)  is  less  than  d2(gj)  a)  >  u 

„c  -  21/V/2P1/2(1-[.)-1/2„t  i  2^  . 


,  where 

c 

(A. 11) 


Hence  we  may,  to  adequate  approximation,  ignore  in  assessing  the 

aliasing  which  occurs  in  the  quadratic  case. 

We  see  from  (A. 5)  that  the  minimum  value  of  N  necessary  to  ensure 
that  A2(oi)  is  leas  than  some  maximum  value  e  is  given  by  N  >  N^,  where 

N2  -  p"1(l-p)e'1  .  (A. 12) 


We  see  from  (A. 3)  that  N2  is  smaller  than  by  the  factor  e.  Hence 
aliasing  is  likely  to  be  less  serious  in  the  quadratic  case  than  it  is  In 
the  linear  case. 

In  order  to  assess  the  implications  of  the  present  model  concerning 
ground-based  observations  of  solar  oscillations,  one  will  need  to  have 
detailed  estimates  of  the  expected  spectrum  (in  particular,  the  spacing 
and  relative  power  of  nearby  lines)  and  the  expected  cloud  cover  at  three 
or  four  observatories  positioned  round  the  world.  It  is  also  desirable 
that  the  present  model  should  be  extended  by  considering  separate  values 
of  p  for  each  of  the  observatories,  and  possibly  by  taking  into  account 
the  correlation  between  cloud  cover  on  consecutive  days. 
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Nevertheless,  we  can  Illustrate  the  results  of  this  model  by  consider¬ 
ing  a  hypothetical  situation.  Suppose  that  four  observatories  are  located 
around  the  world  in  such  a  way  as  to  give  continuous  coverage  (in  the 
absence  of  cloud  cover),  and  that  these  observatories  are  operated  for  25 
days.  Then  N  •  100.  Suppose  that,  for  any  observatory  on  any  day,  there 
is  20%  probability  of  cloud  cover  so  that  p  -  0.80.  We  find  from  (4.1) 
that  (using  v  ■  o)/2tt),  ■  15yRz.  For  frequencies  less  than  this  value, 

(4.5)  shows  that  the  aliasing  amounts  to  0.25%  or  less.  For  frequencies 
above  15uHz,  the  aliasing  drops  off  rapidly. 

Although  it  will  be  necessary  to  make  more  detailed  and  specific 
calculations  to  draw  definite  conclusions,  it  appears  from  the  above  simple 
example  that  it  may  be  possible  to  carry  out  high-quality  studies  of  solar 
oscillations  from  a  chain  of  ground-based  observatories. 

This  work  was  supported  in  part  by  NASA  Grant  NGL  05-020-272,  Office 
of  Naval  Research  Contract  N00014-75-C-0673,  and  the  Max  C.  Flelschmann 
Foundation. 
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